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Already Say About Cosmological Parameters 
in Open and Critical-Density Cold Dark Matter Models? 
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ABSTRACT 

We nse a combination of the most recent cosmic microwave backgronnd (CMB) 
flat-band power measnrements to place constraints on Hnbble’s constant h and the 
total density of the Universe Do in the context of inflation-based cold dark matter 
(CDM) models with no cosmological constant. We nse minimization to explore the 
4-dimensional parameter space having as free parameters, h, Do, the power spectrnm 
slope n and the power spectrnm normalization at = 10. Conditioning on Do = 1 
we obtain h = 0.33 ± 0.08. Allowing Do to be a free parameter rednces the ability of 
the CMB data to constrain h and we obtain 0.26 < h < 0.97 with a best-£t valne at 
h = 0.40. We obtain Do = 0.85 and set a lower limit Do > 0.53. A strong correlation 
between acceptable h and Do valnes leads to a new constraint Doh^'^^ = 0.55 ±0.10. We 
qnote Ax^ = 1 contonrs as error bars, however becanse of nonlinearities of the models, 
these may be only crnde approximations to la conhdence limits. 

A favored open model with Do = 0.3 and h = 0.70 is more than ~ 4cr from the CMB 
data best-£t model and is rejected by goodness-of-£t statistics at the 99% CL. High 
baryonic models (D^h^ ~ 0.026) yield the best CMB x^ fits and are more consistent 
with other cosmological constraints. The best-fit model has n = 0.91?o;o9 and Qio = 
18.01 i; 5 hK- Conditioning on n = 1 we obtain h = 0 . 55 ?o;i 9 , Do = 0.70 with a lower 
limit Do > 0.58 and Qio = 18.0l};5 The amplitnde and position of the dominant 
peak in the best-fit power spectrnm are Apeak = 7Q-7 and ipeak = 260 ? 2 o- 

Unlike the Do = 1 case we considered previously, CMB h results are now consistent 
with the higher values favored by local measurements of h but only if 0.55 ^ Dq ^ 0.85. 
Using an approximate joint likelihood to combine our CMB constraint on Doh^'^^ with 
other cosmological constraints we obtain h = 0.58 ± 0.11 and Do = O.Od^oj®. 
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1. INTRODUCTION 

The ensemble of cosmological data prefers 
best-bet universes which seem to congregate in 
several distinct regions of parameter space (Os- 
triker & Steinhardt 1995, Viana 1996). Among 
the best-bet universes, open models figure promi¬ 
nently and are possibly the favorite candidate 
(Liddle et al. 1996a). This preference is mainly 
due to observational evidence (e.g., Willick et 
al. 1997, Carlberg et al. 1996, Dekel 1997). Fur¬ 
ther motivation for examining Hq < 1 mod¬ 
els is that galaxy cluster baryonic fraction lim¬ 
its seem to be inconsistent with Big Bang nu¬ 
cleosynthesis (BBN) if Ho = 1 and h ^ 0.50 
{h = Ho/WO km s~^Mpc~^). This “baryon 
catastrophe” has led some to believe that Ho < 1- 

Recently, theoretical open universe models have 
been developed. Open-bubble inflation models 
have been developed by Ratra & Peebles (1994), 
Bucher, Goldhaber & Turok (1995), Yamamoto, 
Sasaki & Tanaka (1995). Open hybrid inflation 
has also been considered (Garcia-Bellido & Linde 
1997). 

1.1. What Kind of Open Models We Con¬ 
sider and Why 

CMB measurements have become sensitive 
enough to constrain cosmological parameters in 
restricted classes of models. In Lineweaver et 
al. (1997), (henceforth “paper 1”), we described 
our method and compared CMB data to 
predictions of COBE-normalized critical-density 
universes with Harrison-Zel’dovich (n = 1) power 
spectra. We briefly looked at CDM and flat 
A CDM models by exploring the /i — plane 
and the h — Xq plane. We used predominantly 
goodness-of-fit statistics to locate the regions of 
parameter space preferred by the CMB data. 

In Lineweaver & Barbosa (1998), (henceforth 
“paper 2”), we used a similar technique, again 
in critical-density universes, to explore the 4- 
dimensional parameter space h, Hf,, n and Q. 
We obtained the result that if Ho = 1 (and our 


other assumptions are correct) then the CMB 
data prefer surprisingly low values of the Hubble 
constant; h ^ 0.30. We found that four indepen¬ 
dent cosmological constraints also favored these 
low values in the Hq = 1 models considered. This 
is in contrast to local measurements of h which 
seem to prefer h ~ 0.65 ± 0.15 (Freedman 1998, 
Tammann & Federspiel 1997). 

The Ho = 1 assumption we have made in our 
previous analyses can be considered very restric¬ 
tive since plausible values for Ho in the range 
0.2 ^ Ho 1.0 can change the power spectrum 
significantly. In this work we consider open mod¬ 
els motivated by the question: Does our h ^ 0.30 
result depend on the fact that we limited our¬ 
selves to Ho = 1? Would a favored open model 
{h = 0.7 and Hq = 0.3) be acceptable to the com¬ 
bined CMB data? What pairs of {h, Hq) values 
are compatible with the CMB data? 

There are reasons to believe that Ho < 1 mod¬ 
els will allow higher h values. In paper 2 we found 
that the position of the primary acoustic peak in 
the angular power spectrum is a dominant fea¬ 
ture determining the low value of h. The posi¬ 
tion of the peak is shifted towards higher i values 
in Ho < 1 models and this should have the effect 
of increasing the h values of the best-fit mod¬ 
els. Motivated by this idea and the more general 
idea of increasing the size of the parameter space 
into interesting regions, in this paper we put con¬ 
straints on the cosmological parameters h, Hq, n, 
and the normalization at £ = 10 in the context 
of Ho < 1 CDM models. We assume adiabatic 
initial conditions with no cosmological constant. 
As in paper 1 and 2, we take advantage of the re¬ 
cently available fast Boltzmann code to make the 
parameter-dependent model power spectra (Sel- 
jak and Zaldarriaga 1996). We do not consider 
Ho > 1 models because the code is not yet avail¬ 
able. 

The recent dynamic interplay between theory 
(providing a fast code to make model specific pre¬ 
dictions) and observations (new measurements 
are coming in about once a month) is increas- 
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ing our ability to distinguish models. Major ef¬ 
forts have been and are being put into obtain¬ 
ing flat-band power estimates. The synthesis of 
these efforts is an important step towards a more 
complete picture of the Universe. Since the main 
goal of two new CMB satellites (MAP and Planck 
Surveyor) is to constrain cosmological parame¬ 
ters, it is important and timely to keep track 
of the data’s increasing ability to reject larger 
regions of parameter space and put tighter con¬ 
straints on preferred models. That is the purpose 
of this paper. 

Previous analyses most closely related to this 
work include Ganga et al. (1996), White & Silk 
(1996), White et al. (1996), Hancock et al. (1998), 
Bond & Jaffe (1997), deBernardis et al. (1997). 
Although methods, models and data sets differ, 
in the limited cases where comparison is possible 
we have found no large discrepancies. 

In Section 2 we summarize the method used 
to obtain the results and examine some of the 
special features of open models. In Section 3 we 
present our h — Qq results and in Section 4 we 
compare them to non-CMB results. In Section 5 
we present our results for n, the normalization. 
Apeak and ipeak- In Section 6 we discuss and 
summarize. 

2. METHOD 

2.1. Data and Analysis 

We use a combination of the most recent 
CMB flat-band power measurements to place 
constraints on h, flo, n and the normalization at 
i = 10, Qio (see Section 2.2 for Qiq definition). 
We examine how the constraints on any one of 
these parameters changes as we condition on as 
well as minimize with respect to the other param¬ 
eters. We obtain best-fit values and approximate 
likelihood intervals for these parameters. 

We update the data of paper 2 to include sev¬ 
eral more points: 

• updated Tenerife point (Gutierrez et al. 1997); 
6Teff = 32.5l^°5^ laK at 4 // = 20 


• new MSAM results (Cheng et al. 1997): ST^jf = 
50 9 ^ IJ.K at 4// = 159 and ST^ff = 65j'^jQ /aK 
at 4 // = 263 

• new preliminary CAT results (Baker 1997); 
STgff = 47.3j'jg3 /xK at 4 // = 422 and ST^ff = 
43.2+{D;f at 4 // = 615 

• new preliminary OVRO result (Leitch 1998); 
5Teff = 56l}i ijK at 4// = 537. 

The current CMB flat-band power estimates 
used in this analysis are listed in Table 1 and 
plotted in Figure 1. Since there is much scatter 
in the data, there is much scepticism about the 
ability of the points to prefer any particular re¬ 
gion of parameter space. We showed in papers 
1 and 2 however that a simple analysis of in¬ 
teresting restricted families of models is capable 
of showing substantial preferences for relatively 
small regions of parameter space. The scatter in 
the data is partially deceiving in the sense that 
averaging the data over broader bands in (. re¬ 
duces the scatter and presents a surprisingly co¬ 
herent power spectrum which roughly follows the 
polynomial fit in Figure 1. 

Essentially, we are trying to find the parame¬ 
ters of the model that looks most like the dotted 
line in Figure 1. Figure 2 is an example of some of 
the model power spectra tested. For each point 
in the 4-D parameter space we obtain a value for 
X^(/i, rio, n, Qio). The parameter values at the 
minimum value {Xmin) best-fit parame¬ 

ters. The error bars we quote for each parameter 
are from the maximum and minimum parame¬ 
ter values within the 4-D surface which satisfies 
X^{h,Q,n,Qio) = xLn + I- To display the result 
we project this surface onto the two dimensions 
of our choice. The calculation is described in 
more detail in papers 1 and 2. 

Figure 3 is the first in a series of contour 
plots which illustrates our results. The four con¬ 
tours correspond to Xmin + surfaces where 
= [1,4,9,16]. The interpretation of these 
contours is not straightforward. The conditions 
under which these contours can be projected onto 
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Fig. 1.— Recent CMB observations compared with the best-fit model from Figure 4 The dotted line is 
a sixth order polynomial fit to the data which has a peak amplitude and position: Apeak ~ 77 /xK and 
^peak ~ 260. The grey region represents the ~ Icr contour in Figure 4; that is, the power spectra from 
models within ~ la of the best-fit model are contained within the grey region. The small squares above 
and below the 5 Saskatoon points represent the 7% correlated calibration uncertainty (Leitch 1998). The 
best-fit model has n = 0.91, Qio = 18.0 ^uK and = 0.026. 

an axis yielding 1, 2, 3 and 4cj confidence inter- rors are normally distributed and ii) the model is 

vals is described in Press et al. (1992 p 690) (see linear in the parameters or that a linear approx- 

also Avni 1976). These conditions are: i) the er- imation reasonably represents the models within 
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Fig. 2.— CMB power spectra showing the influence of Qq, ^bh'^ and h. These models are for n = 1, 
Qio = 17 fiK. The dotted line represents the data. It is the same in all panels and is the same as in 
Figure 1. The peak amplitude Apeak depends strongly on h but also on and Hq- The £ value of the 
peak, ipeak, is flo dependent but also mildly h dependent. In the lower right panel Qq = 0.7, h = 0.50, 
= 0.026 fits the data quite well and is very close to the best-fit model for n = 1 models (Figure 3). 


the range of parameters of interest. deBernardis 
et al. (1997) find that the error bars are approx¬ 
imately normal. 

Although the model power spectra are nonlin¬ 
ear in the parameters, an approximation linear 
in the parameters may be able to represent the 
power spectrum near the minima. For exam¬ 
ple, in Figure 2, several families of models are 


plotted. The data have a maximum £eff value 
of 615. The ability of the x^ to discriminate be¬ 
tween models comes almost exclusively from the 
data in the range 2 < £eff ~ 430. The strongest 
nonlinearities in the models (typified by regions 
where the models do not trace each other but 
overlap in complicated ways) are in the £ range 
to which the data is not sensitive. It could be ar- 
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Fig. 3.— Likelihood contours in the h—Qo plane. 
We condition on n = 1 while Qio and are 
free to take on the value that minimizes the 
value at that point. The four contours corre¬ 
spond to Xmin + where = [1,4,9,16] 
(See Section 2.1). The best-ht parameters are 
h = 0.551q;} 9 and flo = 0.70 with Hq > 0.58 (~ 
Ifj). The lightly shaded region represents the age 
constraint 10 < to < 18 Gyr. 


gued that the relevant parts of the models may 
be approximated by the hrst two terms of a Tay¬ 
lor expansion of the power spectrum around the 
best-ht parameter values. The accuracy of this 
linear approximation is a measure of the accu¬ 
racy of the correspondence we would like to es¬ 
tablish between the + 1 contour and the 

Ifj conhdence interval. The ehipticities of the 
contours in Figures 3, 4 and 7 are also a mea¬ 
sure of the accuracy of the linear approximation; 
exactly linear models would give concentric, ex¬ 
act ellipses with identical orientations and posi¬ 
tion angles and with semimajor axes in the ratios 
1:2;3:4. We conclude that the error bars that we 
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Fig. 4.— Same as previous hgure except here 
we no longer condition on n = 1. Recall that 
our error bars are obtained from the projection 
of the Ax^ = 1 contour onto an axis. Thus at 
~ 1(7, h is free to take on any value between 
0.26 and 0.97. At the minimum, h = O.dOj’jQ^^ 
and flo = 0.85 with flo > 0.53. The elongated 
Ax^ = 1 contour means that h and flo are highly 
correlated. This correlation leads to a new con¬ 
straint: 0.45 < < 0.65 which should be 

compared to the constraint on the same quantity 
from cluster baryonic fractions (see Figure 5 and 
Section |4T|). 


derive from the Xmm + 1 contours may be useful 
approximations to Icr conhdence limits, but that 
the Xmm+'^j 9) 16 contours are at best rough 

guides to the 2, 3, and 4 a conhdence intervals. 
Work is in progress to quantify the accuracy of 
the linear approximation. Bond, Jaffe &: Knox 
(1998) have done a preliminary analysis compar¬ 
ing a x^-iainimization analysis of hat-band esti¬ 
mates to a more complete pixel-based treatment. 
Their general conclusion is that our “radical data 
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compression method works...sort of” since the 
minima found by the two techniques agree fairly 
well. 

2.2. Normalization and Definition of Qio 

We normalize in the middle of the COBE 
DMR data {(. = 10) rather than at the edge 
{(. = 2) to reduce the otherwise strong correlation 
between the best-fit slope and normalization. We 
parametrize the normalization at £ = 10 using 
the symbol Qio defined by | 

10(10 + 1)010 = ^11. (1) 

Equation is simply a way to write Cio with 
the added convenience that for an n = 1 pure 
Sachs-Wolfe spectrum {C^ oc + 1))), Qiq 

is equivalent to the power spectrum normalizing 
quadrupole Qrms-PS (see Smoot et al. 1992). 

2.3. Saskatoon Calibration 

We have used the new calibration (Leitch 
1998) for the Saskatoon results whereby the nom¬ 
inal Saskatoon calibration (Netterfield et al. 1995, 
Netterfield et al. 1997) is increased by 5% with 
a correlated calibration uncertainty around this 
new value of 7%. We treat the calibration of the 
5 Saskatoon points as a nuisance parameter 
coming from a Gaussian distribution with a dis¬ 
persion of 7% (rather than the 14% used in paper 
2). In this sense our error bars include an esti¬ 
mate of the Saskatoon calibration uncertainty. 

In general the minimum fits prefer Ugk ~ 
0.86. This can be understood quite easily by ex¬ 
amining Figure 1. The little boxes above and 
below the Saskatoon points are ±7% of the cen¬ 
tral values. Usk = 0.86 corresponds to —14%. 
Moving all 5 Saskatoon points down by ~ 14% 
gives the best agreement with the dotted line, 

®A CMB skymap can be written AT jTo (n) = 
and its power spectrum is Ct = 
X) /{‘it + 1). Note that here is not an ensem¬ 
ble average. 


representing all the data. Thus Ugk ~ 0.86 is the 
preferred value. 

There are 32 data points and in the most gen¬ 
eral case where all 5 parameters (h, flo, n, Qio, 
Usk) are free, there are 27 degrees of freedom (= 
32-5). When both h and Ho are low and thus 
nominally Qq < we set Hf, = ilo, thus cre¬ 
ating purely baryonic models in the lower left 
corners of Figures 3 and 4. Computer limits re¬ 
strict the number of discrete values we can 
test. We have performed all calculations for each 
of 3 values of G {0.010,0.015,0.026}. 

Thus we have explored three 4-D slices of pa¬ 
rameter space. For completeness we have also 
minimized with respect to in the same way 
we have for the other parameters, but the 
minimization is restricted to only three discrete 
values. The results from this highly discretized 
BBN range minimization are indicated by “ * ” 
in the column of Table 2. 

2.4. Physical Effects in Open Models 

Acoustic oscillations of the baryon-photon fluid 
at recombination produce peaks in the CMB 
power spectrum around degree angular scales. It 
is convenient to discuss power spectra in terms of 
the amplitude and the position of the hrst such 
peak; Apeafe and ipeak- For example, the ampli¬ 
tude and position of the polynomial ht to the 
data (dotted line in Figure 1) are Apeafc = 77 /xK 
and ipeak = 260. For the physics of the acoustic 
peaks, see the pioneering work by Hu (1995) and 
Hu & Sugiyama (1995a, 1995b). 

In Figure 2 we plot CMB power spectra to 
display the influence of Oq and h and The 

dotted line is the same in each panel, is the same 
as in Figure 1 and represents the data. These 
models are for n = 1, Qiq = 17 /rK. In Figure 
2 we can see that the peak amplitudes Apeak de¬ 
pend strongly on h and . Higher values of 
lead to larger Doppler peaks due to the en¬ 
hanced compression caused by a larger effective 
mass (more baryons per photon) of the oscillat¬ 
ing fluid. For a given a small h means high 


7 




^peak- Qio variations would raise and lower the 
entire curve while variations in the slope n would 
raise and lower Apeak- 

In Figure 2, when h Ipeak T- Similarly, as 
i, ^peak !• The Oo dependence of ipeak is 
a purely geometric effect. The more open the 
universe, the smaller the angle subtended by the 
same physical size. The main point of Figure 2 
is that lowering h and lowering flo both have the 
same effect of raising ipeak- In paper 1 and 2 we 
maintained that it was predominantly the posi¬ 
tion of the peak that favored low /i in Qq = 1 
models. Here we see that lowering Hq can raise 
ipeak to fit the data, hence h does not have to be 
as low. The behaviour of Apeak and ipeak can be 
tabulated as. 


■^peak T 

hi Ho J, (for Ho ^ 0.6) H^/i^ f ref 

^peak T 

hi Ho J, 


which is to be read Apeak goes up when h or 
Ho go down or when or n go up”. See Hu, 
Sugiyama &: Silk (1997) for more details. 


3. RESULTS FROM THE /i - Ho DIA¬ 
GRAM 

The h — ilo diagram is a convenient framework 
in which to explore and present a combination of 
cosmological parameters. The regions preferred 
by the CMB are shown in Figures 3 and 4. The 
results from these figures are given in the first 
two sections of Table 2 which also contains the 
main results of this paper for h, Hq, n and Qio- 
For each result, the conditions under which it 
was obtained are listed and these conditions are 
relaxed as we move from the top to the bottom. 
Table 2 also lists values and the correspond¬ 
ing probabilities P{x^ <) of obtaining values 
smaller than the values actually obtained, un¬ 
der the assumption that the errors on the data 
points are Gaussian, x^ values and probabilities 


are discussed in Section 5.3. 


3.1. h Results 

For Ho = 1 we get h = 0.33 ± 0.08, which 
is the same low h value we obtained in paper 2. 
The new Saskatoon calibration and the new data 
used here do not change our previous result. 

In Figure 3 we present the likelihood contours 
in the h — Ho plane for re = 1. The minimum is 
at h = O.SSlg'j^g. The minimum x^ value and the 
probability of obtaining a smaller value for that 
model are Xmm = ^1.4 and P{x^ <) = 23.4% 
respectively. Thus the fit is “good”. 

Figure 4 is the same as Figure 3 except we no 
longer condition on re = 1. The best-fit h value 
stays low but higher and lower h values are now 
acceptable at ~ lex. Thus, uncertainty in re plays 
an important role in the inability of CMB data 
to determine h. The banana-shaped Ax^ = 1 
contour can be projected onto an axis to yield an 
approximation to a Icr confidence interval around 
the best-ht value. Thus, 0.26 < h < 0.97 or 
h = 0.40+°;?^. 

A favored open model with Ho = 0.3 and 
h = 0.70 is more than ~ Aa from the CMB data 
best-fit model and can be rejected based on good¬ 
ness of fit at the 99% confidence level. In con¬ 
trast to our previous Hq = 1 results, allowing 
Ho < 1 permits much larger h values and there is 
no longer a disagreement with more direct local 
measurements of h. 

3.2. Ho Results 

Our Ho results are also given in Table 2. In 
Figure 3 (re = 1) we obtain Ho = 0.70 and set a 
lower limit Ho > 0.58. We obtain no upper limit 
because we were unable to test Ho > 1 models. In 
Figure 4 we obtain Ho = 0.85 with a lower limit 
of Ho > 0.53 at ~ Icr and Ho > 0.43 at ~ 2cr. 
Thus the CMB can place important constraints 
on these models. 

If we assume that h ~ 0.65T0.15 (as indicated 
by local h measurements) then the CMB data 
prefer a density in the range 0.54 ^ Ho ~ 0.84 and 
a power spectral slope 0.96 ^ re ^ 1.12. 






3.3. A New Constraint on 


The variations in n and permit the large 
range of h seen in Figure 4. For example, for the 
highest values of /i, n ~ 1.18 and = 0.026 
while for the lowest h values n ~ 0.85 and 
= 0.010. n variations alone are not suffi¬ 
cient to permit very high h values. For example if 
we let n be free but we condition on ^bh'^ = 0.010 
then h remains small: h < 0.42. If we condition 
on = 0.015 then h < 0.51. n and are 
keeping the peak amplitude fit correctly. High h 
values suppress the peak height but this is com¬ 
pensated for by high n and VLbh?- 

The elongated banana-shaped = 1 con¬ 
tour in Figure 4 means that h and Hq are anti¬ 
correlated. In Figure 2 we see that ipeak T when 
h I oi Qo i- Thus high Qq go with low h and 
low Qo go with high h. The = 1 contour in 
Figure 4 traces out this strong anti-correlation. 
To get a constraint on two parameters simultane¬ 
ously we need to look at the Ay^ = 2.3 contour. 
This can be described by = 0.55 ± 0.10. 

This should be compared to the constraint on 
the same quantity from cluster baryonic fractions 


(see Figure 5 and Section |4.l|) . 


We adopt the 2a limits of the empirical fit of Pea¬ 
cock & Dodds (1994) (see also Liddle et al. 1996a) 
and include the n dependence, 

0.222 <r-0.32f--l 
\n 

with the assumption that 0.8 < n < 1.2. 

4.1. X-ray Cluster Baryonic Mass Frac¬ 
tion 

Assuming that clusters are a fair sample of 
the Universe, observations of the X-ray luminos¬ 
ity and the angular size of galaxy clusters can 
be combined to constrain the quantity 
We adopt the range 0.04 < < 0.10 

(White et al. 1993) with a central value of 0.06 
(Evrard 1997). We include the uncertainty in 
the value of Qb by replacing Qb with the BBN 
range [0.010 /i“^, 0.026 h~'^] yielding the limits 
0.10 < Do < 0.65, and the central value 
Do = 0.25. 

4.2. Limits on the Age of the Universe 
from the Oldest Stars in Globnlar 
Clnsters 


^ < 0.293. (3) 


4. NON-CMB CONSTRAINTS IN THE 
/i - Do PLANE 


To view our results within a larger picture, 
we compare them to other cosmological measure¬ 
ments and identify what the CMB constraints 
can add to this picture. The independent non- 
CMB cosmological measurements are summa¬ 
rized below. They are the same constraints used 
in paper 2 (with modifications described below) 
and we now include local measurements of h. 

Peacock &: Dodds (1994) made an empirical fit 
to the matter power spectrum using a shape pa¬ 
rameter F. For Do < 1 models, F can be written 
as (Sugiyama 1995) 


F = /i Do exp 



/ Do T 1 \ 

)_ ■ 


( 2 ) 


Although the determinations of the age of the 
oldest stars in globular clusters are h- and Do- 
independent, they do depend on the distance as¬ 
signed to the globular clusters of our Galaxy. In 
paper 2 we used 11 < to < 18 Gyr with a central 
value of 14 Gyr. We now adopt 10 < to < 18 
Gyr with a central value of 13 Gyr because the 
recent Hipparcos recalibration of the local dis¬ 
tance ladder increases the distance to the glob¬ 
ular clusters (Feast & Gatchpole 1997, Gratton 
et al. 1997, Reid 1997). This lowers the inferred 
ages by about ~ 5 — 10% depending on what val¬ 
ues were used in the calculation of the globular 
cluster distances. 

Age determinations are Do and Aq indepen¬ 
dent but converting them to limits on Hubble’s 
parameter depends on Do and on our Ao = 0 as¬ 
sumption. For a flat Universe with Ao = 0 and 


9 







h 

Fig. 5.— This plot has no CMB information in it. 
The three bands are constraints from three non- 
CMB cosmological measurements discussed in 
Section The “Age” of the oldest stars in globu¬ 
lar clusters: 10 < to < 18 Gyr, the baryonic frac¬ 
tion in “Clusters” of galaxies: 0.10 < < 

0.65 and the matter power spectrum shape pa¬ 
rameter “T”: 0.169 < r < 0.373. The thick con¬ 
tours are approximate Icr and 2a regions from a 
joint likelihood of these three constraints with 
the added constraint from local measurements 
that h = 0.65 ± 0.15. The results: h = 0.60^'^qJq 
and rio = 0.45 ]'^q^ 7. An uncertainty of the bary¬ 
onic fraction 0.010 < Vt^h? < 0.026 has been in¬ 
cluded in the constraints. In flat models the three 
constraints shown favor low values of h ~ 0.40 in¬ 
compatible with local measurements of h. In the 
open models considered here, this disagreement 
disappears. 

to expressed in Gyr, h = (6.52/to). In an open 



h 


Fig. 6.— Approximate la and 2a contours from 
a joint likelihood of the new CMB constraint 
= 0.55 ± 0.10 with the four constraints 
shown in Figure 5. The results: h = 0.58 ± 0.11 
and Oo = O.OOto'is- Popular Hq = 1 and 
(rio = 0.3,/i = 0.70) models are ~ 2a from the 


best fit. 


universe 

2 - Ho 
Ho 

(4) 

where the term in square brackets is 2/3 for Ho = 

1 and goes to 1 as Ho —> 0. The constraint 10 < 
to < 18 Gyr, inserted into Equation provides 
the “Age” constraint on h used in Figures 3, 4 
and 5. 

4.3. Summary of Constraints Used 

The constraints we adopt from cluster bary¬ 
onic fraction, the ages of the oldest stars in glob¬ 
ular clusters, the matter density power spectrum 
shape parameter F and local measurements of h 
are, 


h = 


9.78 

to 


Ho 


1 - Ho 2(1 - Ho)3/2 


cosh 


-1 
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Clusters 

0.10 

< 


< 

0.65 

Age [Gyr] 

10 

< 

to 

< 

18 

F 

0.169 

< 

F 

< 

0.373 

Hubble 

0.50 

< 

h 

< 

0.80 


where the respective central values adopted are 
rio = 0.25, to = 13 Gyr, F = 0.25 and h = 
0.65. 

The first three constraints are illustrated by 
the three bands in Figure 5. The lo" and 2a 
regions from an approximate joint likelihood of 
all four constraints are also shown (see paper 2, 
Section 4.5 for details). The la region yields: 
h = 0.60](qJq and Qo = 0.45](q^ 7. If we con¬ 
sider only the first three constraints the result 
is h = 0.60lg2i and flo = 0.45 ](q^ 7. Thus, 
in rio < 1 models, there is good agreement be¬ 
tween the first three constraints and local mea¬ 
surements of h. This was not the case for the 
flo = 1 universes tested in paper 2 where the first 
three constraints favored lower values; h ~ 0.40 
(notice in Figure 5 that for Hq = 1, /i ~ 0.40 is 
preferred). 

An uncertainty of the baryonic fraction of 
0.010 < Q.bh? < 0.026 has been included in 
both the cluster and F constraints. We have also 
made a figure analogous to Figure 5 but with a 
smaller BBN uncertainty around a higher value, 
0.022 < < 0.026. For this case, the lower 

limits of the “Cluster” and “F” bands are raised, 
thus narrowing the la region. 

4.4. Comparison of CMB and Non-CMB 
Constraints in the h — Q.o Plane 

What does the CMB add to the larger picture 
provided by these non-CMB measurements? 

• Overall consistency; A superposition of Fig¬ 
ures 4 and 5 shows that the ~ la regions of 
CMB and non-CMB overlap for 0.52 < /i < 0.67 
and 0.58 < Qo ^ 0.71. 

• More detailed consistency; The region of 
overlap of the first three constraints and the 
CMB is in agreement with local measurements 
of h. This agreement between CMB, three in¬ 


dependent cosmological measurements and local 
h measurements is non-trivial; in paper 2, al¬ 
though we had agreement between the CMB and 
three independent cosmological measurements, 
the agreement was at /i ~ 0.35 and did not agree 
with local h measurements. 

• New constraint; We find a tight new model- 

dependent constraint = 0.55 ±0.10 which 

favors the higher values of the cluster constraint 
on this same quantity. 

• Preference for high Inside the = 

1 contour (except for a small region to the left 
of the best fit) of Figure 4, the is for 

= 0.026. The consistency between the non- 
CMB constraints and the CMB constraints is 
stronger for higher values of Qbh'^. This improved 
consistency and slightly better fit indicates that 
high is preferred, lending some support to 
Tytler et al. (1997) values. 

• New argument for n ~ 1; Figure 3, where 
n = 1, has a minimum inside the joint likeli¬ 
hood Icj contour of Figure 5. In this sense it 
is more consistent with the combined constraints 
of Figure 5 than are the results of Figure 4. We 
can turn the argument around and say that the 
non-CMB constraints favor n ~ 1 based on this 
consistency. 

• More precise combined constraint: Com¬ 
bining the CMB constraint on Qoh^^'^ with the 
non-CMB contours in Figure 5 we obtain; h = 
0.58 ± 0.11 and Qq = O.OO^^gJg (see Figure 6). 

• Rejection of a favored model; The Qq = 0.3, 
h = 0.70 model is acceptable to the non-CMB 
measurements but is more than ~ 4cr away from 
the best-fit CMB model in Figure 4 and can be 
rejected based on goodness of fit at the 99% CL. 

Fiddle et al. (1996a) have examined open mod¬ 
els with Ao = 0. They consider the shape param¬ 
eter F, bulk flows, the abundance of clusters and 
the abundance of Ly-a systems and the age of 
the Universe. They find a good fit for Uq > 0.35 
and an alarmingly low h value for Uq ~ 1. As¬ 
suming h > 0.6 (as indicated by many recent 
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Qio {[-iK) 



0.6 0.8 1.0 1.2 1.4 


n 

Fig. 7.— Likelihood contours in the n— Qio plane 
with h and Oq as free parameters. The mini¬ 
mum is n = O.DltQQQ, Qiq = 18 . 01}'5 
tice that there is no strong correlation between n 
and Qio as there is between n and Qrms-PS (e-g- 
Lineweaver 1994). 

measurements) they get 0.30 < flo < 0.60. If 
we assume h > 0.6 we obtain 0.53 < flo < 0.75 
from the CMB analysis, and 0.3 < Q.o < 0.70 
from the hrst three non-CMB constraints. Thus 
we hnd consistent but slightly higher allowed in¬ 
tervals for Hq. 

5. RESULTS FOR n, Qio, ^peak AND 

^peak 

5.1. Results for n and Qio 

Our n results are listed in Table 2. Condition¬ 
ing on Oo = 1 we get n = O.OIIquo- Figure 7 
displays our most general result in the n — Qio 
plane and yields n = O.Oltoog- Thus the min¬ 
imum is unchanged and the error bars increase 
slightly in this more general case. The best-ht 


value of n is a robust result in the sense that it 
does not change from the Oq = 1 to the flo < 1 
case. 

For rio = 1 and n = 1 we get Qio = 
17.0^{'o irK. For Hq free and n = 1 we obtain 
Qio = IS.Otj's Conditioning on Oq = 1 we 
obtain Qio = 17.5l}'i And hnally with all 
other parameters free we get Qio = IS.OlJ'g ;uK 
(Figure 7). We can also express this normaliza¬ 
tion in terms of Cio, 

This normalization at the best-ht values for h, 
rio) n and should be compared to the 

slightly higher, more general (but COBE DMR 
only) Bunn & White (1997) normalization which 
is a function of the hrst and second derivatives 
of the power spectrum at £ = 10 . 

5.2. Results for -^peak «5-Od ^peak 

We can get a rough idea of the values of Apeak 
and ipeak preferred by the data from the sixth or¬ 
der polynomial ht shown in Figure 1. This yields 
Apeak ~ 77 and ipeak ~ 260. We can also 
make a more careful model-dependent estimate 
of Apeak and ipeak by looking at the power spec¬ 
trum of the best-ht model in Figure 4 and by ex¬ 
amining the power spectra from models along the 
edge of the = 1 contour. The power spec¬ 
trum of the best-ht model gives us central values 
for Apeak and ipeak while the power spectra of the 
models along the edge of the = 1 contour 
yield error bars on these central values. The re¬ 
sult is Apeak = 76^7 fxK and ipeak = 260 l 2 [|- It 
should be remembered that these results depend 
on the correctness of the family of models we are 
considering {Qq < 1 , Aq = 0 ). 

5.3. and Probabilities P(x^ <) 

The minimum values are given in Table 2 
along with the probability of obtaining smaller 
values under a Gaussian assumption for the er¬ 
rors on the hat-band power estimates. The range 
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of minimum-x^ values and their corresponding 
probabilities are [20.9 < < 24.6] and [18.2 < 

P < 35.2]. These values are “good” and bor¬ 
der on “too good”. The highest x^ values and the 
highest probabilities are when we condition on 
h = 0.50 with Qo = 1 giving substantially worse 
fits than h free. The lowest probabilities are 
when we condition on Oq = 1 and = 0.026. 


We have added the calibration uncertainty in 
quadrature to the statistical error bars on the 
flat-band power estimates. If we were more con¬ 
servative we would add them linearly. In this 
case the x^ values would be even lower and the 
hts even better, i.e., “too good”. 


Figure 2 shows how the CMB power spectra 
vary as a function of In Table 2 we have 

included results for each value of separately. 
The = 0.010 minima have the highest x^ 
values while the = 0.026 minima have the 
lowest and are thus identical to the x^ of th® 
restricted minimization with respect to (see 
Section 2.3). 


rio = 1 models are a subset of the models ex¬ 
amined here. The differences between the IIo = 1 
results reported here and those reported in paper 
2 are small and can be understood by the three 
differences in the analysis. In the present work i) 
we only look at three discrete values of of Qbh'^ (in 
paper 2 we explored the h — Qb plane), ii) we in¬ 
clude 5 more data points, iii) we use a Saskatoon 
calibration 5% higher with a smaller uncertainty. 


6. DISCUSSION AND SUMMARY 
6.1. Review of Results 

We use CMB flat-band power estimates to ob¬ 
tain constraints on h, Hq, n and Qiq in the con¬ 
text of rio < 1 CDM models of the Universe. 
Conditioning on flo = 1 we obtain h = 0.33 ± 
0.08. Allowing Hq to be a free parameter reduces 
the ability of the CMB data to constrain h and 
we obtain 0.26 < h < 0.97 with the minimum at 
h = 0.40. We obtain flo = 0.85 and set a lower 
limit Qo > 0.53. We find a strong correlation be¬ 


tween acceptable h and flo values leading to the 
new CMB constraint = 0.55 ± 0.10. We 

also obtain n = 0.9llQ;Qg and Qio = 18 . 0 lJ ;5 
High baryonic models iVt.bh'^ ~ 0.026) yield the 
best CMB x^ fits which also are more consis¬ 
tent with other cosmological constraints. We find 
that a favored open model with Qq = 0.3 and 
h = 0.70 is more than ~ 4(T from the CMB data 
best-fit model and can be rejected at the 99% CL 
based on goodness of fit. 

6.2. Consistency with Non-CMB Mea¬ 
surements 

In the flat CDM models of paper 2 we found 
that h ~ 0.30. This value was consistent with 
four non-CMB constraints but in disagreement 
with local measurements of h. Considering Hq < 
1 models for ^bh'^ fixed at 0.015 we again find 
h limited to values ^ 0.50 and best-fit values of 
Ho near 1. It is not until we allow = 0.026 
that a much larger interval of h is allowed at the 
~ Icr level. For this most general case, the results 
from the CMB, the same non-CMB constraints 
as used previously and local measurements of h 
are all consistent with h ~ 0.58 and Ho ~ 0.66. 

The CTs that corresponds to our best-fit model 
is £78 = 0.511(q;2o. This is consistent with results 
from independent measurements which favor the 
interval [0.4, 0.9] (Viana Sz Fiddle 1996, Eke et 
ai. 1996). Our Hq results are also broadly con¬ 
sistent with bulk flow measurements which yield 
roughly Hq > 0.4 (Dekel 1997). 

CMB constraints are independent of other cos¬ 
mological measurements and are thus particu¬ 
larly important. The fact that reasonable x^ val¬ 
ues are obtained means that the current CMB 
data are consistent with inflation-based Hq < 1 
CDM models for a broad range of h values. In the 
context of the models considered, the CMB re¬ 
sults are consistent with three other independent 
cosmological measurements and are now also in 
agreement with local measurements of h. This 
consistency was not present in Hq = 1 models. 
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6.3. Review of Assumptions 

The results we have presented here are valid 
under the assumption of inflation-based CDM 
models with Gaussian adiabatic initial conditions 
and with no cosmological constant. We have 
not considered early reionization scenarios or hot 
dark matter. We have also not included any 
gravitational wave contributions which seem to 
make the fits slightly worse without changing 
the location of the best-fit parameters (Liddle 
et al. 1996b, Bond and Jaffe 1997). With only 
scalar perturbations, deviations of the power 
spectrum from power-law behavior is negligible 
in open models (Garcia-Bellido 1997). Super¬ 
curvature modes are not included in the power 
spectra models (Seljak & Zaldarriaga 1996) be¬ 
cause they are in general important only for £ ^ 5 
and thus difficult to measure because of cosmic 
variance. 

It is possible that one or more of our basic as¬ 
sumptions is wrong, or we could simply be look¬ 
ing at too restricted a region of parameter space. 
Topological defects may be the origin of struc¬ 
ture. Using the same data and y^-minimization 
analysis, we find (Durrer et al. 1997) that several 
classes of scalar-component-only global topologi¬ 
cal defect models also produce acceptable fits to 
the data although the goodness of fit of these 
models is not as good as the models we consider 
here. In other words, goodness-of-fit statistics 
from current GMB data have a slight preference 
for the inflation-based models we have consid¬ 
ered over the topological defect models we have 
considered. 

6.4. Future Improvements 

In addition to the h, Hq, n, Qio and Q.}yh? 
considered here, regions of a larger dimensional 
parameter space deserve further investigation in¬ 
cluding Ao, £1hdm, early reionization param¬ 
eters such as Zreion, tensor mode parameters ut 
and T, the inflaton potential r, iso-curvature or 
adiabatic initial conditions and topological defect 


models with their additional parameters. 

The fact that we obtain acceptable values 
in our 4-D parameter space lends some support 
to the idea that we may be close to the right 
model. If the Universe is not well described by 
these models then as the data improve, work like 
this will show poor fits and other regions of 
parameter space will be preferred. 

To increase the parameter-constraining power 
of the measurements, observations need to be 
made in regions of Gspace that have no or few 
measurements. In Figure 1 we can identify these 
regions: 600 < £ < 1200, 20 < £ < 50 and 
180 < £ < 400. More than a dozen on-going 
small-angular-scale experiments continue to fill 
in these gaps (Page 1997) as we await the more 
definitive MAP and Planck satellite results. 

The improvement of non-CMB measurements 
will reduce the size of parameter space we need 
to look at making the model power spectra com¬ 
putations more tractable. For example if 
can be determined to be = 0.024 ± 0.002 
as claimed by Tytler et al. (1997) (or some other 
equally well-constrained value) then a much smaller 
range of the h — Q.i, plane needs to be examined 
and the range of h allowed by the CMB analysis 
will be much narrower. The indeterminacy of n, 
which seems to be measurable solely by the GMB 
and whose error bar has a relatively large con¬ 
tribution from irreducible cosmic variance, will 
remain a dominant factor in the uncertainty of 
GMB parameter estimation. 

We gratefully acknowledge the use of the Boltz¬ 
mann code kindly provided by Uros Seljak and 
Matias Zaldarriaga. We benefited from useful 
discussions with Max Tegmark and Luis Tenorio. 
G.H.L. acknowledges support from the North At¬ 
lantic Treaty Organization under a grant awarded 
in 1996. D.B. is supported by the Praxis XXI 
CIENGIA-BD/2790/93 grant attributed by JNICT, 
Portugal. 
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Table 1: Data“ Used in the Fits and Plotted 
in Figure 1 


Experiment 


reference 


e// 


DMRl 

DMR2 

DMR3 

DMR4 

FIRS 

Tenerife 

SP91 

SP94 

BAM 

Pythl 

Pyth2 

ARGOl 

ARG02 

MAX GUM 

MAX ID 

MAX SH 

MAX HR 

MAX PH 

Ski 

Sk2 

Sk3 

Sk4 

Sk5 

MSAM 

MSAM 

MSAM 

MSAM 

CAT 

CAT 

CAT 

CAT 

OVRO 


Binshaw et al. (1996 

) 

Hinshaw et al. (1996 

) 

Hinshaw et al. (1991 

) 

Hinshaw et al. (1991 

) 

Glanga et al. (1994| 

"lutierrez et al. (1997) 

Glunderson et al. (1995' 

Glunderson et al. (1995' 

Fucker et al. (1997) 


Platt et aJ. (1997) 


Platt et aJ. (1997) 


leBernardis et al. (1994( 

Vlasi et al. (1996) 


Fanaka et aJ. (1996) 


Fanaka et al. (1996) 

Fanaka et al. (1996) 

Fanaka et al. (1996) 

Fanaka et al. (1996) 

>fetterfield et al. (1997) 

Metterfield et aJ. (1997) 

Metterfield et al. (1997) 

Metterfield et al. (1997) 

Metterfield et al. (1997) 

Oheng et aJ. (1996) 


Oheng et aJ. (1997) 

Oheng et aJ. (1996) 

Pheng et ai. (1997) 

Dcott et at. (1996' 



Baker (1997) 


Dcott et at. (1996' 


Baker (1997) 



3 

7 

14 

25 

10 

20 

60 

60 

74 

87 

170 

95 

95 

138 

138 

138 

138 

138 

86 

166 

236 

285 

348 

159 

159 

263 

263 

396 

422 

607 

615 

537 


^ef.f 


±a' 


data 


+6.6 

4.0 
+3.6 
2.8 
+3.4 
3.1 
+ 14.4 
- 1.2 
.+7.r 


27.9 

24.6 

30.8 

1.2 


29.417 7 
32.5+^°5^ 

30.2tf;5 

36.3^6^1® 
f:; + 29.6 
-15.2 
n+15 


55.6^ 


601i3 

+8.7 
8.7 


39.1 


46.8l®2®i 

46.3l?J:« 

49 1+2^® 

^^•-*--16.4 

32.71^°,'^ 

3-L.o_io 9 
49.0+5« 
+7.0 
6.0 

+ 10.0 
- 8.0 
+ 12.0 
- 10.0 
+ 19.0 
-28.0 
+30.5 
-17.0 
+ 13 
-9 

+23.9 
-14.5 
+ 14 
-10 
+ 13.6 
-13.6 
+9.3 
6.3 
+ 19.1 
-13.7 
+ 13.5 
10.1 
+ 14 
-11 


69.0 

85.0 

86.0 

69.0 

40.7 
50 

44.4 

65 

51.8 
47.3 

49.1 

43.2 
56 


(mK) 


“ CMB anisotropy detections reported in publications since 
1994. The Netterfield et al. (1997) points in Figure 1 are 5% 
higher than these numbers due to the Leitch (1998) recalibra¬ 
tion. See Lineweaver et al. (1997) and Lineweaver Sz Barbosa 
(1998) for further details. 
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Table 2: Parameter Results 


Result® 



Conditions^ 


(p(x^ or 


h 

Oo 

n 

Qio(/iK)'^ 


(%) 

Ha = 35 + f 

- 

1 

1 

free 


22.5(24.0) 

33^1 

- 

1 

free 

free 


21.2(22.2) 

45tl3 

- 

free 

1 

17 


22.0(21.9) 


- 

free 

1 

free 

0.010 

22.6(24.7) 

^'^-13 

- 

free 

1 

free 

0.015 

22.2(22.8) 

ce: + 13 
^^-14 

- 

free 

1 

free 

0.026 

21.4(19.5) 

cc + 13 

00_i9 

- 

free 

1 

free 


21.4(23.4) 

30+J^ 

- 

free 

free 

free 

0.010 

21.6(24.3) 


- 

free 

free 

free 

0.015 

21.3(22.8) 

iotll 

- 

free 

free 

free 

0.026 

20.9(21.1) 

H 

+ 1 

O 

- 

free 

free 

free 


20.9(25.4) 

Ua = 0.85lo ig 

free 

- 

1 

17 


22.0(21.9) 

l.OO^Q 28 

free 

- 

1 

free 

0.010 

22.6(24.7) 

0.85{g,2i 

free 

- 

1 

free 

0.015 

22.2(22.8) 

0 70+®-^® 
^'-'-0.12 

free 

- 

1 

free 

0.026 

21.5(19.6) 

0.70f_oi2 

free 

- 

1 

free 


21.5(23.8) 

1-00-0.21 

free 

- 

free 

free 

0.010 

21.6(24.3) 

0.90l Q ;,g 

free 

- 

free 

free 

0.015 

21.3(22.8) 

0.85lg,32 

free 

- 

free 

free 

0.026 

20.9(21.1) 

0.85^0 32 

free 

- 

free 

free 


20.9(25.4) 

~~n— 1 

n — i-i^'^-o.04 

50 

1 

- 

free 


24.6(35.2) 

^•^-^-0.08 

free 

1 

- 

free 

0.010 

21.6(20.1) 

0 Q1+0.10 
'^•^-^-0.06 

free 

1 

- 

free 

0.015 

21.3(18.7) 

0 Ql +0-07 

'j-yf-o.oo 

free 

1 

- 

free 

0.026 

21.2(18.2) 

0 9i+0-i° 
^•^-L_0.10 

free 

1 

- 

free 


21.2(22.2) 

O-Oll^-.Jg 

free 

free 

- 

free 

0.010 

21.6(24.3) 

o-sstVot 

free 

free 

- 

free 

0.015 

21.3(22.8) 

^•^-^-0.09 

free 

free 

- 

free 

0.026 

20.9(21.1) 

0 91+0 29 

free 

free 

- 

free 


20.9(25.4) 

Qto = i7-oi;:S 

free 

1 

1 

- 


22.5(24.0) 

n.btli 

50 

1 

free 

- 


24.6(35.2) 

18.01?-^ 

free 

1 

free 

- 

0.010 

21.6(20.1) 

iT.sj:}® 

free 

1 

free 

- 

0.015 

21.3(18.7) 


free 

1 

free 

- 

0.026 

21.2(18.2) 


free 

1 

free 

- 


21.2(22.2) 

is.otlt 

free 

free 

1 

- 


21.5(23.8) 

18.0tl;° 

free 

free 

free 

- 

0.010 

21.6(24.3) 

18-01o,5 

free 

free 

free 

- 

0.015 

21.3(22.8) 

18.0+1-2 

free 

free 

free 

- 

0.026 

20.9(21.1) 

18.01 + 2 

free 

free 

free 

- 


20.9(25.4) 


“ parameter values at the the minimum values. The results cited in 
the abstract are in bold. The units of Hq are km s~ Mpc~ 

^ “free” means that the parameters were free to take on any values 
within the discretely sampled ranges: 0.15 < h < 1.00, step size: 
0.05, number of steps=18, 0.1 < ^ 1-0, step size: 0.05, number 

of steps=19, 0.49 < < 1.51, step size: 0.03, number of steps=35, 
12.0 < Qio < 25.0 l-tK., step size: 0.5 A^K, number of steps: 26, 
0.010 < < 0.026, only three values: 0.010, 0.015 and|9r926. Thus 

we have examined more than 900,000 models. See Section p.3l for more 
details about conditions. 

^ Probability of obtaining a smaller x^- There are 32 data points and 
the number of degrees of freedom varies between 26 and 28. 

^ QlO ~ QrJrTsT-PS ~ Qflat pure Sachs-Wolfe, n = 1 power spectra 
(see Section |2.2|). The units of Qio are pK 

^ Highly discretized minimization within the BBN range: E 

{0.010,0.015,0.026}. See Section 2.3 for details. 

^ The Ax^ = 1 contour extends to values rio > 1. 
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